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On My Past Writings 


Kit Tyabandha, Ph.D. 
24 October 2002 


Cyberspace 


In 1983 Apple computers had been installed in number at the Ash- 
burton College, New Zealand, where I was doing my Sixth Form as an 
American Field Service exchange student, and Pascal was taught. In 
1993 I worked as a Security Officer for the TelecomAsia Company in 
Thailand. My job was to look after data on the three mainframe com- 
puters in use. It was there that I learnt terms like cyberspace and hackers. 
I became fascinated by this idea, that there is a virtual space tessellated 
into partitions which house a population whose identity is related to 
real people, because to me it is a whole new dimension added to the 
existing ones. Starting doing my master degree in Control and Infor- 
mation Technology at UMIST in 1994, I came in closer touch with the 
cyberspace community through the use of Unix is the first thing taught 
in that course. I voluntarily wrote two articles on the subject, the first 
one is about cyberspace (Tiyapan, 1994) and the second one about its 
security (Tiyapan, 1996d). 


Sociology 


The nature of discoveries and progresses in science is according to 
Bacon (Francis Bacon, 1620) a birth of Time rather than a birth of Wit. This 
is the same idea of percolation and the description he gave is the very 
picture of the theory. According to him major scientific progresses come 
in revolutions which are sparsely distributed in both time and regions. 
There have only been three periods of major progress out of the five and 
twenty centuries over which the memory and learning of men extends, namely 
the Greeks, the Romans and the nations of Western Europe. These are 
narrow limits of time, the periods in between of which are unfavourable 
to development. A discovery or an invention, then, comes as a chance 
accident in the scale of an individual, and as a certainty when looking 
from a distance. 


When the time is right and all the hidden momentum built up, 
theories will come on by itself as a rule. This does not negate the 
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excellence of an individual, but in a society where there are enough 
multitude of individuals the show will always go on, with or without 
a particular genius. This idea can very well explain cases of multiple 
discoveries. According to Kekulé in his Benzolfest speech in 1890, when 
he ascribed his conception of the cyclic nature of Benzene in dreams, 
certain ideas at certain times are in the air and if one man does not enunciate them, 
other will do so soon afterwards. 


To see the relationship of this with percolation it is possible to 
look at two different things in turn, first at the discoveries and then at 
the discoverers. With a unit being that of a discovery the connection to 
percolation is that big discoveries come as connections of other smaller 
and less obvious ones. A theory often has more than one perspec- 
tive, and which one of them comes to the fore first depends much on 
which combination happens to percolate through first. The discoveries 
of Schrédinger and Heisenberg in Quantum Physics can bear witness 
to this both in the combination and the multiple discoveries parts of this 
argument. 


Let us turn our attention now to the scientist and look at the one 
who does the discovering instead. The theory of percolation tells us 
that at the point of discovery he is by no mean the sole integral ingre- 
dient. If he does not do it, then someone else will certainly do. In order 
to see this, I did four simulations for the cell, bond, vertice, and edge 
percolations on a two-dimensional Voronoi network and then another 
four with the same respective blockage of each case but considering 
the inverse phase instead. The number of units considered for the four 
cases are n, = 200, ng = 416, n, = 298, and n, = 426. With the order 
of simulations as described above, at just one step before percolation 
occurs there are respectively 10.6, 13.4, 11.5, 1.1, 11.24, 10.0, 19.4, and 
7.9 percents among the remaining units which will readily trigger the 
onset of percolation. In other words, these are atoms which are able to 
link up existing clusters and form a percolating cluster. 


The formation of mobs is an interesting phenomenon comparable 
with phase change in physics. What happens is that an agglomerate of 
individuals becomes one and a single creature, the underlying mecha- 
nism of which still baffles any effort towards understanding it. Likely 
enough it has got something to do with psychology and the mind. But 
to me at least, the phenomenon is percolative. Having gained some 
acceptance from my previous writings (Tiyapan, 1995, 1996a, 1996b) I 
gave my new work which briefly discusses the mechanism underlying 
the formation of mobs to an editor of the Sakkayaphab journal whom I 
know. At that time a political turmoil unequal in its degree and extent 
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has been going on for five years. Whether by fate or by design there 
has been a successful but tragic use of mob in Bangkok. The word 
mob has joined the list of those synonym to distrust, namely communist 
or, in western community now for that matter, Islam and evil. Whether 
because of this or something else, the article (Tiyapan, 1996c) simply 
and mysteriously got lost; no one would admit having seen it, and the 
translation of another subsequent article of mine (Tiyapan, 1996d) has 
not been without a noticeable negligence. Thus to me distrust is also 
percolative. The list of things one finds over-distrusted without reason- 
able explanation goes on indefinitely, homosexuality, communism, etc. 
The same seems to be the case with bad habits. My father used to teach 
me using the following poem, 


Bad habit gathers by unseen degree 


like brook makes river, river runs to sea. 


Looking back, it could have been the title of that article, on prag- 
matists and idealists, which has somehow convinced the editors into be- 
lieving that it was political which to me is nonsensical. I only meant 
literary, even if at times philosophical. I include it here because it con- 
tains a curve showing a critical emotional transition. 


The formation of the United States, the European Union, or the 
Commonwealth comes from the trust which acts to join countries to- 
gether like glue boxes in Tex. Like all binding forces, trust is mutual and 
spreads in the same way as a growing cluster does. The cluster grows 
bigger as one or more members are added, and it becomes stronger as 
the level of the mutual trust increases. In a similar way, distrust is also 
mutual and also spreads. If I distrust you and you distrust me, I will 
make sure that I remain as far away from you as possible while you 
will certainly avoid me by all means in return. 


Only these two are possible, so there are only two phases to con- 
sider, that of trust and distrust. The relationship where one trusts while 
the other distrusts would not be stable, since the former will soon learn 
to join the latter. Trust forms clusters of one phase, distrust another. 
The size of these clusters vary in a way similar to those in percolation 
of geometrical networks. The strength of the glue is analogous to the 
probability either of becoming or remaining a member of a cluster. 


The rise of dictators, the proliferation of weapons of mass destruc- 
tion, etc, these things I believe are the products of changes of something 
hidden within the underlying structure. Unless we find out what is 
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happening in the background, these things will unavoidably occur. I 
believe that this unseen thing behind the scene is governed by some 
phenomena similar to that of percolation. I think that the key towards 
understanding many unexplainable phenomena is to investigate, in the 
light of the percolation theory, the working of agglomeration of coun- 
tries or states like those of the United States and the European Union. 


The believe that scientific discoveries are a birth of time, rather 
than of wit (cf Larsen, 1993), is the same as the idea of percolation. 
We know because we remember. And all the various discoveries of 
our time together with the knowledge we possess of the past bring us 
closer to another discovery. Scientific discoveries, then, is the collective 
product of humanity rather than property of a single person (cf Merton, 
1965). 


And because all species are also the product of percolation in time, 
our knowledge and consciousness, too, are the product of the universe. 
One may say that it is a personification when we refer to a collective 
noun, for instance a mob, as an individual. But the truth is that, under 
the percolation theory, it is in fact a separate individual without any 
need for the use of a simile. The renormalisation group theory tells 
us that there exists a structure in a bigger scale that behaves like the 
individual components that comprise it. It seems, therefore, that for 
humans these collective beings of ours are still primitive compared 
with each of us as an individual. This is the reason why, whenever 
we come together, we alway make wars. In our case, then, we seem to 
be conscientiously percolated only individually not collectively. In the 
case of bees, on the other hand, it is the other way round. This is why a 
colony of bees does things which make far better sense than a bee does. 
But one can not say that even a colony of bees has consciousness, because 
there seems to be no morals in what it does. I do not know whether 
there are other beings in the universe both the collective and individual 
beings of whom have percolated conscientiously. But I believe they 
exist, in which case they should be more intelligent than us, though 
this is by no means necessarily the case. 


Control Systems 


The following are relevant subjects in control systems. I have one 
first degree and one master degree in the subject. From 1995 to 1999 I 
did a doctorate course in Tokyo but decided to quit after three years. 
Research in this area is still one of the topics that I would like to do in 
the future. One year before leaving Japan my research results must have 
look quite well because on 224 April 1998 my supervisor Professor 
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Katsuhisa Furuta wrote me an email which says, ‘Your reports are very 
nice! Excellent results. I am verymuch [sic] impressed. We can talk the 
results.’, signed ‘Furuta’. And again on 2”4 June 1998 another email, 
‘Dear mr.kit [sic] please come and explain your synchronous motor. 
Where shall we introduce the controller?’, signed ‘Katsuhisa Furuta’. 


Although it is true that control systems is always used in military 
and missiles, one must not forget that it is also important for travelling 
in space and is likely to be of a great help to those space settlers search- 
ing for habitable planets in the future. Control system is used much in 
electrical, mechanical, and chemical engineering, with the typical time 
constant increasing in that order. 


This section is in a way a brief recapture of what I did during my 
PhD study at TIT, Japan. I collected my works in a form of Technical 
Reports and gave a copy each to Professor Furuta who was my super- 
visor then. Sometimes when there was a spare copy left I would give it 
to another senior staff who worked in the Minami 5 building where the 
Furuta lab used to be. The computer files of these reports is no longer 
available even to myself. A great number of figures from simulation 
results in these reports are not reproducible here without the Simulink 
facility on Matlab. The first of such report was dated 16%” July 1998 
but the works it contains started around the beginning of April of the 
same year. I studied systems listed in a book by Khalil (1996). One of 
these systems (Exercise 1.17 (4)) is 


&y = 21 + 22-21 (|21| + |22/), (1) 
bo = —201 + £2 — £2 (\e1| + |22|) +4, (2) 


where u is the input to the system and 2; and 22 the state variables. 
Because the simultaneous equations 2; + x2 — £1 (|%i| + |z2|) = 0 and 
—22%1+2%2—22 (|21| + |e2|) = 0 are ill-formed, there can be no equilibrium 
points. 
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Figure 1 Phase plane of Equations 1 and 2 with no input 


In Figure 1 is shown the phase plane of the system described by 
Equations 1 and 2 when there is no input. The limit cycle in Figure 1 is 
elliptical with the major axis along the x2-axis. When the control input 
is u = sgn(x1; + £2) the phase plane looks like Figure 2 (a) and when 
u = —(@1 + 22) + sgn(a1 + 22) Figure 2 (b). 
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When there is no control Figure 1 shows that the equilibrium point 
at the origin is an unstable focus and there is a limit cycle circling 
around the origin. This limit cycle is slightly larger in the x2 direction 
than in the x; direction. Every trajectory starting from an initial point 
other than the origin goes to and then stays on this limit cycle. With 
the control input u = sgn(x1 +22) there is a discontinuity on the surface 
of the hyperplane s = 21 + #2 = 0. All trajectories still converge to 
the limit cycle although the latter is distorted where it intersects the 
hyperplane. There is no node and every point on the hyperplane and 
inside the limit cycle is an unstable node. All these also hold when 
the control input is u = —(a1 + x2) + sgn(x1 + 22) and there are two 
additional nodes as shown in Figure 2 (b). 


During September and October 1996 I gave a series of seminars 
on a topic related to polytopes of polynomials. The topic I chose was 
recent (Pujara, 1996) and rather difficult for me but for me it was a 
success because the methods introduced was discussed weekly for at 
least a month in a series of the subsequent seminars, and they resulted 
in the idea being successfully applied, in the context of Pulse Width 
Modulation, by one of the students who attended in his Ph.D. work. 
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The following is a recapture of the original seminar I gave. The dis- 
cussions, which was the most interesting part in these seminars, are 
lost. A pseudoboundary is defined to be the set of all polynomials 
in the polytope each of which has at least one zero on the imaginary 
axis. A section of the pseudoboundary corresponding to wo is a poly- 
tope whose vertices lie in the exposed 2-d faces of the given polytope. 
Pujara (1996) in his study of the stability boundary problem gave an 
algorithm to generate all the vertex polynomials of any section of the 
pseudoboundary of a polytope. A polytope is stable if and only if its ex- 
posed edges are stable. Interval polynomials are a hyperrectangle in co- 
efficient space with edges parallel to the coordinate axis. Kharitonov’s 
results is that if a polynomial family consists of interval polynomials, 
then the stability of just four extreme polynomials is both necessary and 
sufficient for the stability of the entire rectangle. Consider the polyno- 
mial f(s,q) = 8s" + ai(q)s"~! + a2(q)s”-? +... + Qn_1(q)8 + Gn(q) which 
produces an r-dimensional polytope P in R"*?. Here a;(q), i =1 to n, 
are real affine coefficients and qg; < q < qj for every q;, 1 <i<r. For 
fixed g, this polynomial is a point in R"*' whose coordinates are the 
coefficients of the polynomials. The pseudoboundary, , is the set of 
polynomials in the polytope P each of which has at least one zero on 
the imaginary axis. The section, 89, of 8 at wo consists of those polyno- 
mials each of which has a zero at jwo. If Z € C"*! is a set of all zeros 
of this polynomial, then wo € W if and only if dg(jwo) € P such that 
g(jwo) = 0. Every B(jwo) is a polytope which has its vertices on the 
exposed 2-d faces of P, which in turn can be explicitly determined. 


Theorem. (Pujara, 1996) The section of the pseudoboundary of a polytope 
P at any frequency wu is a polytope. The vertices of this pseudoboundary lie 
in the exposed 2-d faces of the polytope P. 


The state equations of a dc motor are 


dw 


di 
L— = —-kw — Ri 4 
7 kw — Ri+u (4), 


where i is the armature current, u the voltage, R the resistance, L the 
inductance, J the moment of inertia, and w the angular speed. The 
constant excitation flux ¢ results in the torque ki and the back e.m-f. 
kw. 
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My technical report number 6 (1998) I studied the control of syn- 
chronous machines. In the introduction I wrote to Professor Furuta to 
say that he may use any of the material without acknowledging me. 
I see it fits to describe briefly here some of the results because I ini- 
tiated and carried out the work, as is normally the case with most of 
my researches done while in Japan. The electrical frequency of such 
machines synchronises with the mechanical speed. However it is more 
convenient to express angles in electrical unit, i.e. electrical angle, rather 
than in mechanical unit, i.e. mechanical angle, for the reason that a syn- 
chronous machine can have more than two poles. Whichever is the 
case, one may consider only a single pair of poles and then consider 
the electrical, mechanical, and magnetic conditions, which are associ- 
ated with all the other pole pairs as repetitions of those for the pair 
being considered (cf Fitzgerald et al, 1971). 


The electrical and the mechanical angles of a generator are related 
by the relation 6, = $6, where 6, is the electrical angle of the output, 
Om is the angle of the rotating shaft, and p is the number of poles of the 
machine. The frequency of the voltage is f = $2, Hz, where n is the 
mechanical speed in rpm and % the speed in revolution per second. 
The radian frequency w of the voltage is we = Swm, where wm is the 


mechanical speed in radian per second. 
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Figure 5 Dynamics of the mechanical and the electrical angles 


Figure 5 shows the dynamics of the mechanical and the electrical 
angles of a synchronous machine with three pair of poles. The solid 
line is the mechanical, while the dash-dot line the electrical rotation. 
The mechanical and the electrical rotation of a synchronous machine 
are juxtaposed in Figure 5. In power system engineering one wants 
to control the power and the frequency. Examples of controllers are 
speed governors, actuators, regulators, and signal transducers. Speed 
governors can be described by Af - kK = Tr, where Af is the change in 
frequency, K the gain, and Tz, the load torque. The other controllers 
mentioned are high-gain power-amplifiers which convert things like 
position, oil pressure, and electricity into valve positions. 


The change in load in the power-frequency control of a synchron- 
ous generator can either be a load shedding where the load decreases 
causing the frequency to increase, a generator shedding where the load 
increases causing the frequency to decrease, or a short circuit where the 
generator is suddenly disconnected and as a result the frequency in- 
creases to the maximum value. The effects caused by the load change 
are related to the equation J ce = Tpm —Tpg, where J is the moment of 
inertia of the rotating mass, @ the angular position of the rotating mass, 
Tpm the prime-mover torque, and Tg the electrical torque. 
374 
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In 1998 I presented a paper (Tiyapan, 1998) at the ATAC-98 con- 
ference in Japan. The objectives of that study are to study the effects of 
a sign function input, to study the effects of ¢ in a singularly perturbed 
system, and to study a variable structure singularly perturbed system in 
other words those which have a sign function in the control input. Sin- 
gular perturbation of any system implies that it has ¢; terms multiplying 
the term x and that these terms €; approach zero. One can write down a 
given transfer function Ts) in the state space forms x = E~1Ax+E~'Bu 


5)(s? 2 
and y = Cx +D. For example, T(s) = <2###482.— can be repre- 


sented by the above state equations with 


-8§ -5 1 1.25 0 1 00 0 

4 0 0 0 0 010 0 
ne O oO -11 -2.5]’ Be 2.82843 | ’ 0 0 1 OF’ 

0 0 4 0 0 000 1 


C=[-2.83 1.6 0.71 0.88], and D=0. 
Then if we rewrite E as 


E= 


? 


0 0 
0 0 
1 0 
O € 


coro 


1 
0 
0 
0 


we can study the perturbation effect of the singularly in € as it decreases 
from € < 1 towards e < 1, namely 





e=0.5, 
T(s) = 2(s + 10)(s* + 4s + 29) 
(s? + 8s + 20)(s + 2.2984)(s + 8.7015) 
e=0.1, 
T(s) = 28° + 108s” + 458s + 2900 
(s2 + 8s + 20)(s? + 11s + 100) 
e=0.01, 
T(s) = 28° + 1008s? + 4058s + 29000 
(s? + 8s + 20)(s2 + 11s + 1000) 
é = 0.001, 
T(s) = 28° + 10008s? + 400588 + (2.9 x 10°2.9) 
s* + 19s? + 1010852 + 80220s + (2 x 105) 
< = 0.0001, 
T(s) = 28° + 1000085 + 4000588 + (2.9 x 10°) 


s* +1983 + (1.0011 x 10°)s? + (8.0022 x 10°)s + (2 x 10°)" 
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To summarise the root locus behaviour when « decreases, two complex 
zeros are fixed, that is s* + 4s + 29. Two complex poles described by 
s’ + 8s + 20 remain the same throughout. The zero on the real axis 
moves left. The two poles on the real axis moves toward each other. At 
the point when the moving zero passes the pole on its left hand side 
the root locus changes its characteristics. 


Even as early as November 1998 I have decided that Control Sys- 
tems will not be the only area that I will do researches in, as can be 
seen in the title of the technical report number two (1998) which has 
the words system and control put in brackets. Here among other things I 
considered the system described by the differential equation 


dr daz 
e—~ + — =u 5 
det dt (5) 
. This can be rewritten as the state equations #1 = x2 and e%2 = —a2+u 


which represent a model in the standard form since the second equation 
has a real isolated root when ¢ = 0. By letting the right hand side of 
the state equation be zero, every point on the 2; axis is an equilibrium 
point. A model described by the equations # = f(t,a,z,¢), « € R” and 
6% = g(t,x,z,€), z € R” is said to be in the standard form if and only 
if 0 = g(t,x,z,0) has k isolated real roots z = hA,(t,x), i = 1,2,...,k 
(Khalil, 1996); then these two equations reduces to a quasi-steady state or 
slow model & = f (t,x, h(t, x), 0). 


For the present system the eigenvalues obtained from the Jacobian 


: _ 2a 
matrix A = of 





= i ‘| are 0 and —L. Therefore it has an 
w=(21,0) 0 == E 


equilibrium subspace and the qualitative behaviour of the trajectories 
depends on the values of ¢, that is when « > 0 all trajectories converge 
to the equilibrium subspace, when ¢ = 0 the system degenerates and 
is reduced to a first order system, and when ¢ < 0 then all trajectories 
diverge from the equilibrium subspace. In other words the system is 
stable when € > 0 and unstable when « < 0. By plotting the state 
planes and the time response of the state variables of this system when 
wu = 0, one can see that the system is stable if ¢ > 0, and unstable if 
€ <0. When |e| becomes smaller the response of the system becomes 
more rapid as ¢2 = 4x2 rapidly converges to its root x. = 0 when ¢€ 
approaches zero. The trajectory is a straight line which becomes more 
parallel to the x2-axis the closer € gets to zero. 
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Figure 6 State plane of the system ed?z/dt* + dx/dt = uwithO<e<1 


As € approaches zero and u = 0 the system becomes 22 = 0 or 
1, = 0 which has the solution x1(t) = 21(0) and x2(t) = 0. The singular 
point ¢ = 0 is the point of discontinuity since the solution is stable 
when e - 0+ and unstable when <e > 0°. 


Next study the effect of a feed back with the sign function, u = 
K sgn(C21 + £2), to the system in Equation 5, where C' and K are real 
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constants. When C = K = 1, simulation results showed that the re- 
sponse of the controlled system changes more rapidly as |e| becomes 
smaller. In other words, the response becomes faster as € approaches 
zero. This increase in speed of the response as € decreases in magnitude 
applies for the unstable case where € < 0 as well as when é€ > 0. 


; 





je 





Figure 7 State plane of ed?x/dt? + dx/dt = sgn(x1 + 22) whene > 1 


The boundary layer or the transient period of the system remains 
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the same either with or without the input, that is t of approximately 
[0,0.05] and [0,0.4] for « = 0.01 and 0.1 respectively. For ¢ > 1 this 
boundary layer extends beyond one second. With the feed back input, 
the state plane trajectory reaches the line x2 = 0 during the transient 
period and then gently slides along it. 


The state equations with the input become #; = a2 and e%2 = 
—x2 + sgn(a1 + a2). When ¢ = 0 the latter equation above becomes 
0 = —a%2 + sgn(x1 + x2) which gives the root £2 = sgn(# + %) that can 
not be written in the form %2 = A(t,Z). This means that the equation 
has no isolated root and consequently one can not isolate the fast mode 
from the slow mode by introducing a new variable y = 22 —h(t,21). But 
for the purpose of finding the boundary layer model, let h(t,21,22) = 
sgn(x; + a2). Then y = x2 — h(t,a1,%2) = v2 — sgn(z1 + #2) and as a 
consequence #2 = y + sgn(z1 + 22). Introduce a new time variable + 

dy dr 


: _ dy 1 oa : = 
obtained from ¢9; = ge Or GF = = When 7 = 0, that is T = 7) + 





ye tdt = —*. This new time variable 7 is generally known as the 
stretched time variable. Then the equation ¢ [s# + 4 sen(a1 + 22)| =-y 
becomes we + sgn(x1 +22) = —y. Substituting x2 = y +sgn(a1 + 22) 
gives su =-y- es sgn(ai+y+sgn(21+22)), and therefore the boundary 
layer model is se =-y- — sgn(x1 +22 +1) when 21 + z2 > 0, is se = 
-y- 4 sgn(x1 +22) when 21+ 22 = 0, or is su =-y- sgn(x1 +22 —1) 
when 2, + 42 < 0. When £, + 2 > 0 then #2 = +1 and the reduced 


problem becomes #; = +1. Likewise when %; + %2 = 0 then % = 0 and 
“, = 0, and when #1 + 2 < 0 then #2 = —1 and “; = —1. 


The discontinuous nature of the signum function makes the anal- 
ysis of the variable structure control system difficult. This is due 
to the fact that one can not analytically find Osan ate) | oognel) where 
i =1,2,...,n and z(t) is an n-tuple vector. One way to overcome or go 
around this problem is to use a numerical approximation for the sign 


function where necessary. 


Let the input be u = —(#1 + #2) + sgn(#1 + 22), then simulations 
show that the system has two timescales, which is the characteristic 
of a singularly perturbed system. Having a fast response within the 
boundary layer and a slow response elsewhere, the boundary layer 
becomes narrower and the response faster as € approaches zero from 
above (€ + 0+). Also, numerical studies shows that this boundary layer, 
in seconds, is approximately [0,0.05], [0,0.23], and [0,0.54] respectively 
when ¢ is 0.01, 0.1, and 1. With this input, the equations of the system 
become #, = x and eg) = —a — 242 + sgn(x1 +22). The last equation 
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becomes = 0 — 2; — 222 + sgn(x1 + #2) when ¢ = 0, in other words 
t= —52) + 5sgn(2x1 +22). To find the boundary layer, change the 
variable x2 to y = 42 — 22|c=0 = X2 + $21 - 5sgn(x1 +22). Therefore rz. = 
y — 42, +35 Sgn(a1 +22) =y—4a1+4s8gn (x1 +y — 521 + $sgn(z1 + 22)). 
Which is essentially that 2. equals y — 521 + }sgn(y + 521 + 4) when 
v1 +02 > 0, y—$a1+5 sgn(yt+$21) when 21+22 = 0, and y—$a1+5 sgn(y+ 
eal - 5) when #1 + #2 < 0. Substitute this into the equation for e%2 


above to get the boundary layer model € [3 7 5oH + Zsgn(a1 + 2)| = 


—2y ore oy Sy ae $2} — ¢sgn(x1 +22) + $sgn(21 +2£2)| = —2y. Use the 
stretched time variable 7 introduced above. Then su = —fy - fa - 


+sgn(x1 +22), which is the same as saying that su equals —3y — 4a; —4 
when 2; + 22 > 0, —3y — fa1 when 2; +22 = 0, and —3y — $x; + ¢ when 


a1 +2%2 <0. 


Now for a variable structure control, let the hyperplane be de- 
scribed by s = cx; + £2 whose derivative is $ = c#; + ad2. Let the 
input be u = —K sgns = —K sgn(cx1 + cx2) for K € R*. Choose a Lya- 
punov function as V = 38°. The state equations become #; = #2 and 
£2 = —2%2 — +(ct1 + 22). The derivative of the Lyapunov function is 


1 K 
then V = 88 = (cx; +22) (cx, — ta. — £ (cx + 22)) = Caja. - SH + cx3- 


Bs — E41 sen (cx; + x2) — *#2(cx; +22). In order to observe the effect of ¢ 
on the stability of the controlled system, let c = 2.7 and K = 3.3. From 
plots of V against x; and x2 one may see that whether V be positive or 
negative depends upon the value < takes. For example, there exists x 
which gives V > 0 when ¢ is 1.7 and 0.08, while if ¢ = 0.3 then V < 0 


for all possible values of x. 





The original system (Equation 5) has the eigenvalues at \; = 0 


and Ay = —12.5, with the corresponding eigenvectors at v; = | and 
v2 = one | respectively. The system is linear since A is a con- 


stant matrix and does not depend upon z. Also A is not a stability 
matrix since the condition Re; < 0 is not satisfied. The equilibrium 
point at the origin is not unique as it has already been pointed out 
earlier that every point on the a ,-axis is an equilibrium point. The 
x,-axis is a nontrivial null space of the matrix A. The change of 


1 -0.0797 

: — y-l = = 

variables z = M~°a, where M = [v1 v2] = E 0.9968 cE leads to 

Mz=AMz= Y and therefore A = eC : The solution 
0 -12.5]’ O -12.5]" 
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to this system is 21(t) = z1(0) and z2(t) = z2(0)e~ 12" = z0(0)e*2". That 
A is no stability matrix can be seen by trying to solve the Lyapunov 
equation PA + A’ P = —Q, where Q is a positive definite symmetric 
matrix. The Lyapunov equation for this system yields no solution. 


From simulation results, when € < 0 the controlled system is un- 
stable, while for ¢ > 0 increasing K will widen the range of € in which 
the controlled system is stable. That is, the larger the value of the gain 
kK of the control input, the more robust with respect to € this control 
scheme. 


Note also that one way to stabilise the origin of the system is to 
solve the equation A? PE+E? PA—(E? PB+S)R'(B? PE+S")+Q =0, 
or equivalently F’PE + E’ PF —-E’PBR'B'PE+Q-SR''S" =0, 
where F = A—BR™'S™, Q and Rare symmetric positive definite matri- 
ces, and the system is written in the form E¢ = Ax + Bu. For example, 


: 12.62 1 
letting Q = R= E=Tand $= 0 leads to P =| 1 ayiasa he 


feedback input is then u = —Gz, where G = R-1(B’ PE+S"). The gain 
matrix willbe G=[1 0.1194]. Then the closed loop eigenvalues can be 


computed from (A— BG)V = EVT, where V=[v1. w],T = ke 7 I 
2 


i; is the i eigenvalue, and v; its corresponding eigenvector. In solving 
the Riccati equation the Frobenius norm of the relative residual matrix 
is 3.14 x 10°16. Simulations on Simulink show that when « = 0 the 
system is unstable. There exists a domain of attraction outside which 
the trajectory does not converge. 
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